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Abstract. With a direct demodulation method, we have 
reanalyzed the data from the COMPTEL observation of 
PKS0528+134 during the 1993 March flare in 7-rays. Our 
results show that during the flare 7-rays were detected at a 
level approximately 2.4-3.8 times higher than the observed 
intensity in two earlier COMPTEL observations VP and 
VP 1 in the energy range 3 MeV to 30 MeV. The 3-30 
MeV time variability of the flux follows well the trend as 
observed by EGRET at higher energies. No convincing 
excess can be found around the position of PKS0528+134 
in the energy range 0.75 MeV to 3 MeV, which indicates 
a spectral break around 3 MeV. The detections and non- 
detections in the four standard COMPTEL energy bands 
are consistent with the earlier reports given by CoUmar et 
al., while the feature that 7-rays of the quasar still kept 
on flaring at energies down to 3 MeV is clearly found. 
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1. Introduction 

The quasar PKS0528-I-134 was first detected at 7-ray en- 
ergies by EGRET aboard the Compton Gamma-Ray Ob- 
servatory ( CGRO ) (Hunter et al. 1993). Before this, it 
had been known as a bright radio source with a flat ra- 
dio spectrum, a optical source with a mean brightness of 
m„=20 and little polarization (Fugmann & Meisenheimer 
1988) and also an X-ray source (Bregman et al. 1985). It 
has a redshift of 2=2.07 (Hunter at al. 1993) and, along 
with the strong 7-ray sources Crab and Geminga, is lo- 
cated near the Galactic anticenter. Following the detec- 
tion of PKS0528-H134 by EGRET, the COMPTEL data 
were analyzed and PKS0528-I-134 was discovered at soft 
7-rays (0.75-30 MeV) by CoUmar et al. (1993a). CoUmar 
et al. (1993b) also reported the preliminary result from 
COMPTEL viewing periods and 1 that spectral breaks 
were found at energies around 10 MeV when comparing 
with the spectra detected by EGRET. PKS0528-f 134 was 
seen to flare during viewing period (VP) 213 by EGRET 



and 7-rays (> 100 MeV) were detected at a level approx- 
imately three times higher than the observed intensity in 
earlier observations VP and VP 1 (Mukherjee et al. 
1996). The quasar was also detected by EGRET in its 
'low' states during the other observations include viewing 
periods 2.1, 36.0 , 36.5, 39, 221, 310, 321, and 337 (Nolan 
et al. 1993). Comparisons of the results from COMPTEL 
VP 213 with those from EGRET were carried out by CoU- 
mar et al. (1997) in their analyses of PKS0528-I-134 using 
the first 3.5 yr of COMPTEL observational data (VP to 
VP 337) . They subdivided the data into two parts accord- 
ing to the luminosity of the source: viewing periods 0, 1 
and 213 as high states, and all the other observations as 
low states. They found that during the low states the com- 
bined spectra can be fitted by a single power-law represen- 
tation, while a spectral break at energies around 10 MeV 
exists for the combined spectra of the high states. How- 
ever, their results for the fiux show that PKS0528-I-134 
had high luminosity in the 3-10 MeV band and, at the 
upper COMPTEL energy band (10-30 MeV), the quasar 
only had normal luminosity during VP 213. This feature 
is odd and stimulates us to reanalyze the flaring quasar. 

Considering the poor statistics of the data from 
COMPTEL VP 213 which covered only 6 days, we adopt 
the direct demodulation method which owns the capabil- 
ity of reconstructing objects from incomplete and noisy 
data to make a spatial analysis upon the flaring quasar. 
We introduce the instrument and the analysis method in 
Sect. 2 and give our results in Sect. 3. The results show 
clearly the detection of PKS0528+134 in the energy range 
3 MeV to 30 MeV. We also compare the intensity ratios 
and the light curves with those derived from EGRET and 
give our flaring spectrum. In Sect. 4, we give a discussion 
and draw the conclusion that the quasar kept on flaring 
in the energy range 3 MeV to 30 MeV during COMPTEL 
VP 213. 

2. Instrument and analysis method 

COMPTEL is one of the four instruments aboard CGRO 
and is sensitive to 7-rays in the energy range 0.75 MeV 
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to 30 McV. COMPTEL consists of two detector arrays. 
An incident 7-ray photon is first Compton scattered in 
the upper detector and then interacts with the lower de- 
tector. Two orthogonal coordinates (x,'0) which describe 
the direction of the scattered 7-ray can be obtained from 
the positions of interaction occurring in the two detectors 
and the Compton scattering angle (p in the upper detec- 
tor can be calculated from the measured energy deposits 
of the photon in the two detectors. The incident direc- 
tion of the photon is then known to lie on a projected 
circle (with its center at (x-l^) and with radius called 
'event circle') on the sky for the ideal case of total absorp- 
tion in the lower detector. The scattering direction (x,'0) 
and the Compton scattering angle (p constitute a three- 
dimensional data space in which the spatial response of the 
instrument is cone-shaped. For the details about COMP- 
TEL see Schonfelder et al. (1993). 

Inversion methods are needed for source reconstruction 
in three-dimensional data space and a direct demodulation 
method is chosen by us. The direct demodulation method 
has already been successfully applied to imaging analysis 
of scan observation data of slat collimator telescopes, e.g. 
the scan imaging of CygX-1 by the balloon-borne hard X- 
ray coUimated telescope HAPI-4 (Lu Z.G. et al. 1995) and 
reanalyzing the EXOSAT-ME galactic plane survey (Lu 
F.J. et al. 1996). It can also be applied to analyze obser- 
vational data from a coded-mask aperture telescope (Li 
1995), rotating modulation telescope (Chen et al. 1998) 
and Compton telescope (Zhang et al. 1997). 

The principle of direct demodulation (Li & Wu 1994) is 
to perform a deconvolution from the following correlation 
equation under proper physical constraints: 

P'! = c (1) 

where P' = P^P, c = P^d, P is the point spread func- 
tion matrix of COMPTEL and P^ its transpose, d is the 
observational data and / the intensity distribution of the 
unknown sky and background. The physical constraints 
can be the upper and lower limits of intensity. The pro- 
cess of direct demodulation consists of two main steps. 
The diffuse background is derived first by iteratively solv- 
ing the correlation equation under constraints of continu- 
ity. Then, the intensity distribution of the object sky is 
obtained by solving the correlation equation again under 
constraints of the produced diffuse background as a lower 
limit. After subtracting the diffuse background, the true 
intensity distribution of the sources in the object sky and 
the non- uniform component of the background are left. In 
this work iteration calculations are carried out by using 
the Gauss-Seidel algorithm and error estimates derived 
by the bootstrap technique. 

3. Results 

^From the beginning of the all-sky survey of COMPTEL 
PKS0528-I-134 has been in the field of view several times. 



During observation VP 213 (March 23 - 29, 1993) the 
source was located about 9.0° off the COMPTEL point- 
ing direction. With the direct demodulation method, we 
have analyzed the data from VP 213 in the four standard 
COMPTEL energy bands. 

Images of PKS0528-I-134 were obtained only in the 
3-10 MeV and 10-30 MeV bands (see Fig. la and Fig. 
lb). In the 0.75-1 MeV and 1-3 MeV bands, no con- 
vincing excess around the position of PKS0528-I-134 was 
found. Fig. la and Fig. lb show clearly the existence of 
a source at the position of PKS0528-I-134. In addition to 
the quasar, the Crab nebula and pulsar Crab also exist 
in the two skymaps and these two 7-ray sources arc re- 
solved exactly. The flux density values of PKS0528-I-134 
are found to be (16.6±5.95)xl0-5 ph cm'^ s'^ in the 
3-10 MeV band and (7.71±3.03)xl0-'^ ph cm'^ s"! in 
the 10-30 MeV band. The relative statistical significance 
is about 2.79cr and 2.54c7. The flux values of the Crab in 
the two skymaps are (3.57±0.68)xl0~'' ph cm^^ s^^ in 
Fig. la and (9.57±2.57)xl0-5 ph cm'^ s'^ in Fig. lb. 
Therefore, although the quasar became the brightest 7- 
ray source during its flare at energies above 100 MeV, it 
was still weaker than the Crab at MeV energies. The flux 
values of PKS0528-M34 obtained by Collmar et al. (1997) 
arc (15.2±4.7)xl0-^ ph cm^^ g-i ^^^^ ^yj^^v band 
and (3.3±1.7)xl0-5 ph cm^^ g-^ in the 10-30 MeV band. 
Our 3-10 MeV flux is close to their result and 10-30 MeV 
flux is about 2.3 times higher. Both the 10-30 McV fluxes, 
nevertheless, are consistent with each other within the er- 
ror bars. Collmar et al. (1997) also gave an upper limit of 
the 0.75-1 MeV flux and a very marginal 1-3 MeV flux, 
which indicates that PKS0528-I-134 was very weak and 
could hardly be detected at energies below 3 MeV. Our 
non-detection of PKS0528+134 in the 0.75-1 McV and 1-3 
MeV bands confirm this indication. The other structures 
in both the skymaps reflect the non-uniform background 
from instrument and the statistical fluctuations of the ob- 
servational data. 

We list in Table 1 the fluxes detected by EGRET and 
by COMPTEL at energies above 3 MeV during three high 
states of PKS0528-I-134. Fluxes marked (*) are our results 
and fluxes of PKS0528-I-134 for COMPTEL VP and VP 
1 are from Collmar et al. (1997). All the EGRET results 
in Table 1 are from Mukherjee et al. (1996). With these 
fluxes, we have calculated the intensity ratios for every 
two high states observed by EGRET and by COMPTEL 
respectively, and the results are listed in Table 2. We can 
see from Table 2 that the 10-30 MeV intensity ratios are 
close to those detected by EGRET at energies above 100 
MeV and the intensity ratios of VP 213/VP 1 and VP 
213/VP in the 3-10 MeV band are also high, which sug- 
gest that the flaring of PKS0528+134 was in phase in both 
the energy intervals 3-30 MeV and > 100 MeV. The non- 
detections in the lower energy bands show its cessation. 
The COMPTEL 3-10 MeV and 10-30 MeV light curves to- 
gether with the EGRET one (Mukherjee et al. 1996) are 
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shown in Fig.2. In Fig.2 our fluxes of PKS0528+134 for VP 
213 are used and all other data points in the COMPTEL 
3-10 MeV and 10-30 MeV light curves are from Collmar 
et al. (1997). From Fig.2 one can sec that the COMPTEL 
3-10 MeV and 10-30 MeV energy bands follow well the 
intensity trend as observed by EGRET above 100 MeV. 

The COMPTEL spectral results together with the si- 
multaneously measured EGRET spectra of PKS0528-hl34 
for VP 1 and VP 213 are shown in Fig. 3a and Fig. 3b. 
The spectrum of the quasar in Fig. 3a is from Collmar ct 
al. (1993b). The spectrum at energies above 30 MeV in 
Fig. 3b is from Mukherjee et al. (1996) and the other two 
data points in this figure arc our results. The COMPTEL 
spectral results shown in Fig. 3b are consistent with a sim- 
ple extrapolation of the measured EGRET spectrum. The 
spectral break in the 10-30 MeV band for COMPTEL VP 
1 (see Fig. 3a) does not appear in the flaring spectrum for 
VP 213 (see Fig. 3b) even down to the energy of 3 MeV. 
The non-detection of PKS0528-I-134 at energies below 3 
MeV implies that a spectral break may occur around 3 
MeV in Fig. 3b. Such a spectral break is also required by 
the results of Collmar et al. (1997). 

4. Discussion and conclusion 

We have carried out an imaging analysis upon the quasar 
PKS0528+134 during its flare with the direct demodu- 
lation method. All our results agree with those given by 
Collmar et al. (1997) except for a difference in the 10-30 
MeV flux. The 3-10 MeV flux obtained by us is close to 
that by Collmar ct al. (1997) and our 10-30 McV flux is 
about a factor of 2.3 higher. Taking into account, however, 
the significance of the detections and the errors of the de- 
rived fluxes, one can see that the two values for the 10-30 
MeV flux are in agreement. The flux difference may result 
from the different treatment of the background when an- 
alyzing the data. In our analysis, we search for the diffuse 
background only from the data and use them as lower limit 
when solving again the modulation equation. The non- 
uniform instrument background and the intensity distri- 
bution of the real sources are left in the final map. In this 
way, we need no longer consider the background prior to 
estimating intensities of the sources because both of them 
can be derived from the data simultaneously. We note 
that there is some structure near PKS0528-I-134 in Fig. lb 
and no other structures aroimd PKS0528+134 in Fig. la, 
which imply that the background around PKS0528-I-134 is 
relatively clean in the 3-10 MeV band but dirty and noisy 
in the 10-30 MeV band. The complicated background in 
Fig. lb may therefore lead to the difference in the 10-30 
MeV flux when using different analysis methods. 

Ahhough PKS0528+134 has been fully studied by 
Collmar et al. (1997), some important features for the flar- 
ing quasar are found by us in our reanalyzing COMPTEL 
data of VP 213 with the direct demodulation method. 
The flaring quasar PKS0528-I-134 was detected with high 



luminosity in the energy range 3 McV to 30 McV and 
not detected in the lower COMPTEL energy bands; The 
intensity ratios and the light curves measured by COMP- 
TEL at energies above 3 McV arc consistent with those by 
EGRET; A spectral break at energies around 3 MeV is re- 
quired from the combined spectrum measured simultane- 
ously by COMPTEL and EGRET during VP 213. These 
features suggest that PKS0528-I-134 still kept on flaring 
in the 3-30 MeV band at similar level as measured by 
EGRET above 100 MeV and no flaring below 3 MeV. The 
inverse Comptonization (IC) models, which can be subdi- 
vided into the self generated synchrotron Compton model 
(SSC model; e.g. Maraschi ct al. 1992, Bloom & Marscher 
1993) and the external radiation Compton model (ERG 
model; e.g. Dermer & Schlickeiser 1993, Sikora et al. 1994, 
Blandford & Levinson 1995) according to the origin of low- 
energy photons scattered to 7-energies by the relativistic 
jet electrons, are generally used to explain the production 
of high energy 7-rays from quasars. To explain the spectral 
break for the flaring spectrum of PKS0528-I-134, however, 
we need a combined ERG / SSC model (Bottcher & Coll- 
mar 1998) which introduces a high-energy electron pair 
population with a low-energy cutoff in its Lorentz-factor 
distribution (Bottcher et al. 1997). Our result for the flar- 
ing spectrum is consistent with this model. Our detections 
and non-detections of the flare at MeV energies put hard 
constraints on the relative low-energy cutoff, thus helping 
to further understand the emission process that produces 
the flare. 
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Table 1. Fluxes of PKS0528+134 measured by COMPTEL and EGRET in three high states. 



VP 


EGRET 


COMPTEL 




> 100 McV 


10-30 McV 3-10 McV 




(xlO-^'ph/cm-^s-^) 


(xlO-'^ph/cm-^s-^) (xlO-^ph/cm-^s-i) 





12.9±0.9 


3.2±1.0 4.4±2.5 


1 


8.5±0.8 


3.0±1.0 < 7.2 


213 


30.8±3.5 


7.71±3.03* 16.60±5.95* 



Table 2. Intensity ratios of PKS0528+134 measured by COMPTEL and EGRET. 



Intensity ratio 


EGRET 


COMPTEL 




> 100 Mo\' 


10-:!0 Me\' 


3-10 AloV 


VP 1/VP 


0.66±0.08 


0.94±0.43 


< 1.64 


VP 213/VP 


2.39±0.32 


2.41±1.21 


3.77±2.53 


VP 213/VP 1 


3.62±0.54 


2.57±1.32 


> 2.31 



Skymap of VP213 (3-10MeV) 
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Fig. 1. a The direct demodulation map of VP 213 in the 3-10 MeV band, b The direct demodulation map of VP 213 

in the 10-30 McV band. Apart from Crab and PKS0528+134, the other source features emerged in these two maps 
are generally insignificant and most of their statistical significances are below 2a. 
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Fig. 2. Time history of PKS0528+134 as measured by COMPTEL in the 3-10 MeV and 10-30 MeV bands and 
by EGRET above 100 MeV (Mukherjee et al. 1996). Data points with square symbol are our results and all other 
COMPTEL data points are from CoUmar et al. (1997). The error bars are la and the upper limits are 2a. 
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Fig. 3. a Energy spectrum of PKS0528+134 measured by COMPTEL and EGRET during VP 1. b Energy spectrum 
of PKS0528+134 measured by COMPTEL and EGRET during VP 213. 



